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Abstract. We explore the complementarity between two tomographic probes 
of the universe: baryon acoustic oscillations (in the galaxy power spectrum) 
and weak gravitational lensing. The galaxy power spectrum characterises the 
density fluctuations, whereas the weak lensing shear power spectrum is a direct 
measure of the potential fluctuations. We find that photometric measurements of 
baryon oscillations alone do not provide very tight constraints on the dark energy 
equation of state parameters, partially due to our uncertain knowledge of the 
galaxy clustering bias. Weak lensing, on the other hand, is adversely impacted by 
the uncertainties of the probability distribution of photometric redshift errors. A 
joint analysis of the two, however, is more robust to these uncertainties and leads 
to a remarkable improvement over the results of either probe alone. 

Forecasts of cosmological constraints with baryon oscillations and weak 
lensing are provided for four proposed multiband imaging surveys in combination 
with measurements of the cosmic microwave background from Planck. In 
particular, we find that the joint analysis of galaxy and shear power spectra with 
the Large Synoptic Survey Telescope can tighten the la error bounds on the dark 
energy equation of state (at the pivot expansion factor a p = 0.63) and its rate of 
change, respectively, to 0.016 and 0.16 (marginalized over 131 other parameters). 
With supcrnovae and cluster counting as well as higher-order statistics of the 
same galaxy and shear data, one can further improve the constraints. 

To allow greater flexibility in utilising our forecasts, we will make our 
forecasting tool available publicly. 

Keywords: cosmological constant experiments, gravitational lensing, power 
spectrum, surveys galaxies 



1. Introduction 

The apparent accelerated expansion of the universe revealed by type la supernova 
(SN) distances 0, H 

suggests the existence of an unknown component - dark 
energy - of the universe that drives the acceleration with its negative equation of state 
(EOS). Observations of the large-scale structure and cosmic microwave background 
(CMB) also point to the scenario that the universe is nearly flat with more than 70% 
of its content in dark energy H, IE 0j El E| ■ 

Since it is not directly observable, one has to infer dark energy properties 
from its impact on the distance-redshift and growth-redshift relationship. Aside 
from SNe, baryon acoustic oscillations (BAOs) in the galaxy power spectrum and 
weak lensing (WL) can be used to measure the angular diameter distance. The 
BAO technique relies on the standard ruler of the sound horizon at the last 
scattering surface 0, 0, i which has been accurately measured by the Wilkinson 
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Microwave Anisotropy Probe ( WMAP) [j| and will be further improved by future CMB 
experiments. WL measures the distance from its geometric lensing kernel as well as 
the shape of shear pow er sp ectra, and it is sensitive to the growth-redshift relationship 
ElQIIlIIlIIlIIlIIllilllll. Cluster counting is another useful technique that 
exp loit the volume-redshift relationship and evolution of cluster number density [e.g. 
l22l |23j , although it requires a very good understanding of the distribution of the mass- 
observable relation including its scatter [3 an d the dispersion of the cluster mass 
function. 

The prospects for measuring the dark energy EOS parameters with BAO and WL 
have been studied in great detailjlljl^^mll BEllIllIJilllllllllllli 

In the absence of systematic errors 
and uncertainties of photometric rcdshift (photo- z) error distribution, a deep and wide 
multiband WL survey can achieve statistical precisions of ~ 0.1 on the dark energy 
EOS parameters wo and w„, , where the EOS is parametrized as w — wq + w a (l — a) 
and a is the scale factor |45j. Shear systematics and uncertainties of photo- z error 
distribution can severely degrade WL constraints 0, 0, 0, 0] . In addition, the 
lack of complete understanding of the matter power spectrum on small scales due to 
nonlinearity and baryonic effects [Hfl l^ll Ir51 IHij makes it more difficult to fully 
harness the statistical power of the shear power spectrum on those scales pifij . 

Given a deep photo-z survey, BAO constraints on wq and w a are generally weaker 
than those of WL, because the rms photo-z error <j z filters out radial information 
on small scales and because we do not have an accurate description of the galaxy 
clustering bias 0, 0, US EH E-1 • When only the auto power spectra are considered, 
BAO results are sensitive to the uncertainties of the photo-z error distribution [ifllliij . 
However, the cross power spectra between redshift bins, as we discuss below, can self- 
calibrate the photo-z bias Sz and rms error a z . As such, tomographic measurements 
of BAOs are much less sensitive to the uncertainties of the photo-z error distribution. 

Although we treat dark energy phenomen olog ically with a loose connection to the 
cosmological constant or quintessence 0, E3, Ea j one should bear in mind that 
dark energy could also be a manifestation of our incomplete knowledge of gravity. 
For instance, gravity in an extra dimension can accelerate the cosmic expansion 
Matching different theories to the same distance-redshift relationship does not 
guarantee a match of the growth redshift relationship 20, 62, 63, 64, 65, 66] . Hence, 
it is important to employ multiple techniques, especially WL, to not only strengthen 
constraints on dark energy properties but also potentially distinguish between different 
theories. 

We focus on incorporating the two emerging techniques BAO and WL with 
CMB measurements. Both BAOs and cosmic shear have been recently detected and 
used to constrain cosmological parameters ElElESESElEilliEi More 
encouragingly, it has been demonstrated that one can achieve an rms photo-z error 
a z < 0.04(1 + z) with deep photometries and sufficient spectroscopic calibrations |75|. 
With a homogeneous population of luminous red galaxies (LRGs), the rms photo-z 
error can be even smaller, and the error distribution is fitted well by two Gaussians 
centred at nearly the same true redshift [j^ . This has led to very recent detections of 
BAO signatures from photo-z LRG samples in the Sloan Digital Sky Survey (SDSS) 

and MegaZ-LRG catalogue [?8| at comparable statistical significance levels as 
the detections from spectroscopic samples in SDSS [T^ and the 2 degree field galaxy 
rcdshift survey 

To combine BAO with WL in a single survey, one has to properly account for 
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the galaxy-shear correlation [79j- We further develop the work in [33| to include the 
uncertainties of the photo- z error distribution and simple additive systematics. Photo- 
z errors, even if their distribution is known perfectly degrade cosmological constraints 
by suppressing the signal. For WL, the lensing kernel itself is much broader than the 
rms photo-z error <r z , so that it is not overly sensitive to photo-z errors. However, 
the uncertainties of the photo-z error distribution introduce uncertainties to the shear 
signal. Given that future WL surveys can measure shear power spectra to percent 
level at I ~ 1000, WL results will be sensitive to the uncertainties of the photo-z error 
distribution. For BAO, the situation is the opposite. Its kernel is defined by the photo- 
z error distribution (assuming that the widths of the tomographic bins match the rms 
photo-z errors), so that BAO power spectra and their cosmological constraints are 
more sensitive to photo-z errors. But the distinctive dependency of the tomographic 
galaxy power spectra on photo-z errors is actually an advantage that can be taken to 
constrain the photo-z error distribution. Consequently, BAO results are less prone to 
the uncertainties in the photo-z error distribution. 

Photo-z errors can have long tails and sometimes can be catastrophically wrong 
due to the confusion of spectral features (especially the Balmer break for galaxies 
at z < 1 and the Lyman break for galaxies at z > 2) or the lack of features in 
the wavebands [for a review, seeE3|. Such catastrophic errors can be mitigated by 
incorporating galaxy apparent magnitude and size priors or by supplementing near- 
infrared data |3, El EH EH • The outliers have been report edly limited to less than 
10% or even completely eliminated in some cases [TH . l82l Isil l85| . 

We idealise the photo-z errors to be Gaussian, so that they can be parametrized 
with the rms a z and bias Sz. The Gaussian simplification allows for a pedagogical 
understanding of the effect of the uncertainties in the photo-z error distribution. 
The BAO constraints on the photo-z error distribution will be weakened when more 
parameters are used to model the distribution j8||. However, since WL cannot self- 
calibrate the photo-z parameters, even very loose error bounds on these parameters 
from BAO will keep WL constraints on dark energy from being degraded without 
bound. Spectroscopic calibrations are crucial for both properly modelling the photo-z 
errors and constraining their distribution, but, it will be very challenging to fairly 
sample galaxies for spectroscopy over the entire area and depth of future imaging 
surveys. Therefore, alternative means of constraining the photo-z error distribution, 
such as photo-z galaxy power spectra and cross-correlations between spectroscopic 
and photometric galaxy samples E3|j wm be valuable to these surveys. 

Photo-z surveys can be used in many ways to probe the universe, and a 
single technique can utilise a variety of statistics. We specifically investigate the 
complementarity between tomographic galaxy and shear power spectra. Higher- 
order statistics can be very useful [l^, \mL |42( , but to properly combine them with 
other statistics, one needs to carefully account for the correlations between different 
statistics. 

The rest of this paper is organised as follows. We define the tomographic galaxy 
and shear power spectra as well as the CMB power spectra in section [3 and discuss 
the Fisher information matrix for error estimation in section |3 Our forecasts on 
cosmological constraints are given in section 0] for four proposed multiband imaging 
surveys (as far as BAO and WL are concerned): Dark Energy Survey! (DES), Large 

| Sec http:// www.darkenergysurvey.org/ 
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Synoptic Survey Telescope§ (LSST), Panoramic Survey Telescope and Rapid Response 
System|| (PS4), and Supernova/ Acceleration Probed (SNAP). We discuss the results 
and conclude in section |SJ 



2. Observables 



Our observables consist of angular maps of galaxy number density n(0), shear 7(0), 
and CMB temperature and polarisation fluctuations. The tomographic galaxy and 
shear power spectra are the covariance of galaxy number density and shear in Fourier 
space, as are the CMB temperature and polarisation power spectra to temperature and 
polarisation. In what follows, we loosely refer to these power spectra as observables 
as well, despite that using power spectra as true observables in the forecast requires a 
slightly different formalism. 



2.1. Galaxy and shear power spectra 

With the Limber approximation |88l Isflj . the angular power spectrum of n(6) and 
7(0) can be written as [33j | 

^2 



pXY 

ij 



9ir 2 f 

M = ^3 J dzH(z)D A (z)W t x (z)W[(z)A 2 s (k;z), (1) 

where lower case subscripts correspond to the tomographic bins, upper case 
superscripts label the observables, i.e., X = g for galaxies or 7 for shear, H(z) is 
the Hubble parameter, Da(z) is the comoving angular diameter distance, A^(fc; z) is 
the dimensionless power spectrum of the density field, and k = £/D A (z). BAO and 
WL do not necessarily use the same binning. In other words, the bin number is defined 
for each technique separately. The window function is 

b(z)^— X = g 

n -' l:l = 4 3 n m H$p A (z) r°° mWD^z') (2) 

2 H(z) ac J z Z D A {z>) 7 ' 

where b(z) is the linear galaxy clustering bias, and Q m and Hq are, respectively, the 
matter fraction at z — and Hubble constant. The galaxy redshift distribution rii(z) 
in the ith tomographic bin is an average of the underlying three-dimensional galaxy 
distribution over angles, and the mean surface density fii is the total number of galaxies 
per steradian in bin i. The distribution rii(z) is broader than the nominal width of 
the tomographic bin (defined in photo- z space) because of photo- z errors. 

We only include galaxy power spectra on largely linear scales, e.g., the scales 
of BAOs, in our analysis, so that we can map the matter power spectrum t o galaxy 
power spectrum with a scale- independent but time-evolving linear galaxy bias |90ll91| . 
One may extend the analysis to smaller scales with a halo model to describe the scale 
dependency of the galaxy bias and, in fact, can still constrain the scale-dependent 
galaxy bias to 1% level [33j. 

Cast in the form of (|TJ and (J2J, it seems as though the WL kernel W?(z) is 
sensitive to the matter fraction Q m (or the matter density oj m — £l m h 2 , where h is 



§ See http://www.lsst.org/ 

|| See http://pan-starrs.ifa.hawaii.edu/ 

% See http://www.snap.lbl.gov/ As its name suggests, SNAP is also designed to measure SN redshifts 
with spectroscopy. 
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the Hubble constant in units of 100 kms -1 Mpc -1 ). Potential (or spatial curvature) 
perturbations, for which inflation often predicts a nearly scale-invariant behaviour, 
are more fundamental than density perturbations. WL shear and CMB temperature 
fluctuations are also driven directly by the fluctuations in the gravitational potential 
rather than the density. Therefore, we recast the angular power spectra in terms of 
the potential power spectrum: 

with 



P i W = / dzHD A WfWfA^(k), (3) 



2a b 



X 



cHD A J z Hi D A (z') 

We have dropped the argument z in © and 10} where there is no confusion. One 
may choose to normalise to either the density, e.g., the rms density fluctuation as 
within a radius of 8ft-~ 1 Mpc at z = 0, or the potential, e.g., the CMB temperature 
power spectra, but one should propagate the factors of o; m through to ensure that 
gravitational lcnsing responds to potential directly. Moreover, precise measurements 
of CMB power spectra have made it more desirable as well as convenient to normalise 
to the potential fluctuations. 

Observationally, the power spectra will have contributions from galaxy shot 
(shape) noise n~ l {^ ms n~ l ) and systematic noise N* s : 

P^(i) - P* Y {1) + S$ Y (< + , (5) 

where S^ Y an d are Kronecker delta functions. For galaxies, g rms = 1, and, for the 
shear, 7 rms ~ 0.2 is due to the intrinsic shape of galaxies and measurement errors. 
The systematic noise 7V S ^ S may arise from (the residuals after correcting for) the 
spatially varying dust extinction, pho tometr y offsets, point spread function (PSF), 
instrumentation effects, and so on |92l I93I l94j. We have assumed in (JSJ that N^L is 
additive, uncorrelated between BAO and WL, and independent of tomographic bins 
or multipoles. This is a great simplification over the additive errors considered in [4flj ]. 
which shows that the degradation to WL results by the systematics is generally not 
bound. The systematic noise we include is among a sub-set of (unrealistic) additive 
errors whose degradation effect can be shown mathematically to be less than ~ 100% 
[iif . As such, we have effectively made an aggressive assumption that all other types 
of systematics are well under control. 

Previous forecasts for BAO use onl y t he three or two-dimensional galaxy auto 
power spectrum within each redshift binjllliililliilMIIlElIi- This is a good 
approximation for spectroscopic surveys, because the intrinsic correlation (due to the 
large-scale structure) between two redshift bins is small if the bins are wide enough. 
For photo-z surveys, the cross-bin power spectra Pf^{tj are no longer negligible due to 
photo-z errors. The left panel of figure H shows that P^(£) (broken lines), depending 
on the binning and photo- z errors, can have comparable amplitudes to that of the auto 
power spectrum Pf^(t) (solid line). Not including the covariance between observables 
in different bins, i.e., the cross power spectra, generally leads to more optimistic error 
estimates. 
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Figure 1. Left panel: Angular galaxy power spectra Pf^{i) with bin i centred 
at z p = 1.07 and bin j centred at z p = 0.72 (dotted line), 0.88 (dashed 
line), 1.07 (solid line), 1.26 (dash-dotted line), and 1.48 (long-dash-dotted line). 
The shaded area indicates la statistical errors of the auto power spectrum 
AP ifM = i^ + 0-5)f sky ]~^ 2 Pif(e), where / sky is the sky coverage, for LSST. 
The right end of each power spectrum curve is set by an f max for each tomographic 
bin to reduce contaminations from nonlinearity. Right panel: Sensitivity of the 
observables fiifijPf^(£) to the photo-z bias <5z. The derivatives are evaluated 
at I = 40 to avoid contributions from baryon oscillations that are prominent at 
I > 100. Though not shown, nifijPf^{i) is also sensitive to the photo-z rms 
error. The grey vertical bands correspond to the relevant tomographic bins in 
photo-z space, while the upper ticks mark the positions of photo-z parameters in 
true redshift space. For demonstration purpose, we use 15 tomographic bins from 
z p = 0.15 to 3.5 in photo-z space and 17 photo-z parameters uniformly distributed 
from z = to 4 in true redshift space. The widths of the bins are proportional to 

(1 + *p). 



Aside from statistical necessity, tomographic galaxy power spectra will also be 
very useful for parameter estimation with photo-z BAO experiments, because they 
provide valuable information about photo-z errors. The reason is that, with Limber 
approximation, the cross power spectra Pf^{£) are determined by the overlap of galaxy 
number density distributions rii and rij in true redshift space and thus sensitive to 
photo-z bias 5z and rms error a z . This is demonstrated in the right panel of figure ^ 
where derivatives of the new observables + nifijPf^(t) with respective to photo-z bias 
parameters (see section [3] for details) are given at £ = 40. The extents of the relevant 
tomographic bins are indicated by grey vertical bands, and the locations of the photo-z 
parameters are marked by upper ticks. The photo-z bias at an arbitrary redshift is 
linearly interpolated between two photo-z bias parameters. 

We refer to the galaxy bin centred at z p = 1.07 (the subscript p following z 
signifies photo- zs) as bin i in the following discussion. An increase of the photo-z bias 
Sz at z = 0.82 means that the redshifts of galaxies around z = 0.82 are more likely 
to be overestimated. As a result, more of these galaxies will be assigned to bin i. 
With no change to the contributions of galaxies from other redshifts, the auto power 
spectrum nfPjf 6 (solid line) increases correspondingly. The product of the fiducial 

+ It is argued in l48l that nifLjP?? ' (£) is the fundamental observable for WL. Our forecasts for 
cosmological parameters are not affected by this change of observables, but the constraints on photo- 
z parameters are slightly improved. 
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galaxy clustering bias and linear growth function evolves less than 1% in this redshift 
range, so that the change in the mixing of galaxy number fluctuations at different 
redshift is sub-dominant to the effect of increased fii . Similarly, an increase of Sz at 
z = 1.29 leads to fewer z ~ 1.29 galaxies assigned to bin i and thus lowers nfPf^ . 

The behaviour of the auto power spectrum nfP^f as a function of the photo- z 
bias Sz near the nominal centre of bin i depends on the overall galaxy distribution 
n(z). An increase of Sz at z — 1.06 (by linear interpolation, it increases the photo- z 
bias in the range 0.82 < z < 1.29) will shift some galaxies from bin i — 1 to bin i and 
some from bin i to bin i + 1. Since the fiducial n(z) peaks at z = 1, the net result is 
more galaxies in bin i and, hence, a positive response of nfPj^ to Sz at z = 1.06. 

The cross power spectra are more sensitive to the tail of the photo-z error 
distribution, because they are determined by the overlap between two tomographic 
bins in true redshift space. Although the increase of Sz at z = 0.82 leads to fewer 
galaxies in lower-redshift bins (dotted line and dashed line), the increased overlap 
between these bins and bin i, i.e., bin i getting more z ~ 0.82 galaxies, results in a 
boost to nifijPf^ for j < i. The increase of Sz at z = 1.06 reduces the number of 
z ~ 1.06 galaxies assigned to the lower-redshift bins, so that UifijPff decreases. The 
behaviour of the cross power spectra between higher redshift bins and bin % (dash- 
dotted line and long-dash-dotted line) can be explained similarly. 

For our forecasts, we use CMBFAST version 4.5.1 j9f| to calculate the matter 
transfer function at z — 0, and then apply the linear growth function and Peacock & 
Dodds fitting formula (PD96) to obtain the nonlinear matter power spectrum 
at any redshift. Details are given in |Appcndix A| and |Appcndix B|_ To reduce 
contaminations of the small-scale nonlinearity and barvonic effects |32L I37I IsiL Iftol 
we set a smallest angular scale in our analysis. Separately, to avoid 
possible effects of dark energy clustering that have not been considered in our analysis, 
we also discard galaxy and shear statistics on very large angular scales [l9|. This 
limits the multipoles 40 < £ < 2000 for WL. The BAO kernel is rather narrow, so it 
is more appropriate to set the criterion in k space. We require that the dimensionless 
power spectrum Ag(k;z) < 0.4 in each tomographic bin and that 40 < t < 3000 for 
BAO. The maximum multipole for a galaxy bin centred at z p follows approximately 
£ max = 340z p + 346^2 at 0.15 < z p < 0.9, 446 - 658z p + 908z2 at 0.9 < z p < 2.1, and 
3000 at z p > 2.1. Setting the minimum multipole to 100 or reducing the maximum 
multipole by 1000 for both BAO and WL does not severely impact our forecasts (see 
section EJ). 

2.2. Cosmic microwave background 

CMB temperature and polarisation power spectra are calculated with CMBFAST for 
scalar modes: 

Cf Y = 2(2tt) 4 J dlnfeA|(fe)Af (k)Aj (k), (6) 

where X — T for temperature or E for E-polarisation, and Af(k) is the radiation 
transfer function |97| . With noise, the power spectra become 

cr = cr + S XY «,;.y„.;.X,) , (7) 

where #t>,i is the full width half-maximum of the beam in channel i, X rms ^ is the noise 
per pixel in temperature or polarisation measurements, and Ba is the beam window 
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function [97j,|98j. We model the CMB experiment after Planck with instrumentational 
characteristics tabulated in We assume a sky fraction of f s ^ y — 0.8 after 

a foreground cut, and limit the multipoles 2 < I < 2000 for temperature and 
10 < £ < 2500 for polarisation. 

CMB temperature and polarisation are nearly uncorrelated with the galaxy 
distribution e xcept for the late-time integrated Sachs- Wolfe (ISW) effect on large 
angular scales |l0Cj . Since we do not use BAO or WL at £ < 40, it is straightforward 
to combine CMB with BAO and WL. 



3. Error estimation 



Assuming that the likelihood function of the true observables, e.g., the map data, is 
approximately a multivariate Gaussian around the maximum, one can use the Fisher 
information matrix to estimate the lower bound of errors of the parameters that could 
be inferred [To l Hoi Ho3 Hoif . 



3.1. Fisher matrix 

The Fisher matrix of parameters {p a } is given by [10 

^ = \^C-^C-^+ d -f C -^, (8) 

2 dp a dp/3 apa opp 

where C and \x are, respectively, the covariance and mean of the true observables. 
For our purpose, the observables have zero mean, and the contributions to the Fisher 
matrix from different spherical harmonic modes are separable, so that 

F a/3 -/ 8ky ^^-TrC, —C t (9) 

with {C e )* Y = P* Y {£) for galaxies and shear (BAO+WL), and {C t ) uv = Cf v 
for CMB. The factor / s k y approximates the effect of incomplete sky coverage with 
a reduction of modes. This is valid for a contiguous survey with roughly the same 
angular size in all directions. The total Fisher matrix is a direct sum of the BAO+WL 
and CMB Fisher matrices. 

Equation @ can be reduced to 

F - - fe j>+i>7s£«r , f£ do) 

where Qi is a column vector containing unique elements of Ci, and Mf equals the 
covariance of Qi in the Gaussian case |33| . From IjlOfl , it may seem possible to discard 
the cross power spectra in error forecasts. However, the Fisher matrix (JHJ and (|10|l 
are based on the multivariate Gaussian likelihood of the true observables, e.g., 

- 21n£(r^) = In \Cf\ + nJiCf^m + const, (11) 

where ng is a column vector of the galaxy angular distribution in multipole space, 
whereas the likelihood becomes 

- 21n£'(n,) = £ [lnC^ + r^(C^)- 1 ] + const, (12) 

i 

if the cross power spectra are discarded. The two likelihood functions C and C can 
differ considerably, when, for instance, photo-z errors induce strong correlations of 
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the galaxy number density fluctuations between the bins. Discarding the correlations 
between the obscrvables may lead to underestimation of the errors on cosmological 
parameters [e.g., EH. However, in the case of BAO, the galaxy cross power spectra 
are crucial for calibrating the photo-z error distribution, which is, in turn, critical 
to measuring the dark energy EOS parameters accurately with the WL technique. 
Consequently, the joint BAO and WL constraints on dark energy can be weakened if 
one does not utilise the galaxy cross power spectra (see the footnote of table OJ • 

To use the power spectra as true observables, one can replace /x and C in (jHJ) with 
the power spectra and their covariance, respectively. In doing so, one also assumes 
that the likelihood function of the power spectra is a multivariate Gaussian, which is 
generally not the same as assuming a Gaussian likelihood function for the map data. 

3.2. Parameters 

Our parameter set includes 11 cosmological parameters, 40 galaxy bias parameters, 40 
photo-z bias parameters, 40 photo-z rms parameters, one systematic noise parameter 
for BAO, and one systematic noise parameter for WL. There are 133 parameters in 
total for the joint BAO and WL analysis. When the BAO and WL techniques are 
considered separately, their respective irrelevant parameters are held fixed. 

The cosmological parameters are the dark energy EOS parameters Wq and w a , the 
matter density L0 m , the baryon density Wb, the angular size of the sound horizon at the 
last scattering surface 9 S , the equivalent matter fraction of curvature f2j<, the optical 
depth t to scattering by electrons in the reionized intergalactic medium, the primordial 
helium mass fraction Y p , the spectral index n s of the primordial scalar perturbation 
power spectrum, the running of the spectral index a, and the normalisation of the 
primordial curvature power spectrum A%(ot ^A^, in matter era) at k — 0.05 Mpc -1 . 
We adopt the three- year WMAP results |9j for fiducial values of the parameters: (u>o, 
w a , u m , wt, e s , fi K , r, F p , n s , a, A^) = (-1, 0, 0.127, 0.0223, 0.596°, 0, 0.09, 0.24, 
0.951, 0, 2. Ox 10~ 9 ). The reduced Hubble constant h — 0.73 and the present equivalent 
matter fraction of dark energy f2x = 0.76 are implicit in our paramctrization. 

The photo- z bias and rms error parameters are implemented according to 48]; 
we evenly space the parameters from z = to 4 and linearly interpolate the values 
at arbitrary redshifts. Note that the photo-z parameters are assigned in true-redshift 
space and that they are not attached to any galaxy or shear bins. The galaxy clustering 
bias is parametrized in the same way. This scheme is based on the expectation 
that photo- zs and the galaxy bias behave more or less smoothly from one redshift 
to another. It is found that the WL forecasts for a DES-like survey become rou ghly 
convergent when the redshift interval of the photo-z parameters is less than 0.15 |48| . 
The convergent redshift interval may be different for other surveys, but it will be very 
expensive to achieve the same precision of the photo-z calibration via spectroscopy at a 
much finer interval. We expect that ambitious photo-z calibrations will be carried out 
at ~ 0.1 redshift intervals (LSST photo-z calibration white paper*) and that similar 
measurements of the galaxy clustering bias will be available. This amounts to 40 
parameters each for the galaxy bias, photo-z bias, and photo-z rms over < z < 4. 
The actual calibration requirement may be less stringent if one combines BAO with 
WL (see section^J. The photo-z and galaxy clustering bias parameters extend beyond 
the nominal cut-off (photometric) redshift of the surveys to accommodate the galaxy 
distribution of the last tomographic bin in true redshift space. Our forecasts are 

* See [http://www.lsst.org/Scieiiee/Phot-z-plan.pdfl 
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Table 1. Priors and additive systematics a . 



Case 


<T P (<^) b 


op (In b) 


Ni ys 




Opt. 


0.05 <t z 


0.15 


lO" 8 


10-8 


Pes. 


0.2 cr z 


0.3 


2 x 10~ 7 


2 x 10~ 7 



a The additive systematic errors are defined in 151 , 
b (rp(a z ) = V2a P (6z). 



not affected by this extension, because the low galaxy number density and poor 
measurements reduce the statis tical significance of the data around the cut-off redshift. 

We assume b — 1 + 0.84 z 106], Sz = 0, and a z oc (1 + z) for the fiducial model. 
The galaxy distribution m in each tomographic bin is determined from the overall 
galaxy redshift distribution n(z) and photo-z parameters 

m =n(z)P(z® l ,z® l ;z) = n{z)I(z% ti , z® ti ; z)/I(0, oo; z), (13) 

where z^ i and Zp i define the extent of bin i, and P(Zp i: z^ ^ z) is the probability of 
assigning a galaxy that is at true redshift z to bin i. The unnormalised probability 
/(a, b; z) is given by 0] 

(z p -z- 5zf 



I(a, b; z) = — I dz p exp 

v27rer. 



2<r? 



(14) 



and the normalisation 7(0, oo,z) is enforced because photo- zs are non-negative. Note 
that even though the probability distribution of photo- zs at a given true redshift is 
assumed Gaussian, the reverse is not true. In other words, the Gaussian assumption 
is flexible enough to allow for modelling of more complex galaxy distributions in 
tomographic bins |48(. 

The derivatives in l|10(l are taken two-sided numerically with steps in parameters 
A(™ o ,uv,w m ,w b ,0 s ,ft K ,T,Yp,n s ,a) = ±HT 3 x (20,40,6,1,3,10,10,20,5,5). The 
derivative with respect to A|. is analytic. The steps of the galaxy clustering bias 
parameters are ±0.005&, and those of the photo-z parameters are ±0.005(1 + 2;). When 
differentiating with respect to the photo-z parameters, we keep the overall underlying 
distribution of galaxies, i.e., n(z), invariant, because the true redshift of a galaxy 
cannot be affected by photo-z algorithms. This approach is complementary to that of 
[iol I44I lifll Isfj ] , in which the number of galaxies in each tomographic bin is fixed or 
tightly constrained while n(z) is allowed to vary. The two approaches are found to 
be consistent for WL |49|. For BAO, n(z) will be largely degenerate with the galaxy 
clustering bias. Given that our BAO results on the galaxy bias are roughly consistent 
with those in 40], which achieve ^2% constraints for a z < 1.3 and / s k y = 0.1 survey 
with ap(6z) = 0.01(1 + z), the two approaches can be consistent for BAO as well. 

We treat the additive systematics of BAO and WL as nuisance parameters. Since 
they can be determined very well from the data, the additive systematics affect the 
constraints mostly by boosting the covariance |49(. Our forecasts will not improve 
even if N& ys and A^ s are fixed. This is somewhat misleading, because even a very 
slight modification to the form of the systematics can lead to unrestrained degradation 
to cosmological constraints for WL |4jj]. Therefore, one must adequately calibrate and 
control all sorts of systematics in order to achieve the constraints we forecast. 

We incorporate independent measurements of galaxy bias and photo-z parameters 
as Gaussian priors in the Fisher matrix. Since the results depend on these priors and 
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Table 2. Survey specifications. 





Area 




Depth 


n 




0"z 


Survey 


deg 2 


Filters 


r mag (10 <r) 


arcmin - 2 


z* 


(1 + *) 


DES 


5000 


griz 


24.1 


8 


0.3 


0.07 


LSST 


20000 


ugrizy 


26.5 


50 


0.5 


0.05 


PS4 


20000 


grizy 


24.1 


8 


0.3 


0.07 


SNAP 


1000 


9 bands a 


26.6 


100 


0.7 


0.04 



a From optical to near-infrared. 



the systematics to some degree, we list our assumptions in table^for both pessimistic 
and optimistic cases. Note that we use a subscript P for priors, i.e., o~p(p a ), to 
distinguish them from constraints o-(p a ). For Gaussian photo-2 errors, a prior of 
crp(6z) = 0.05o\j translates to 400 fairly sampled spectra for calibration around the 
redshift of the photo-2 bias parameter. Current d eterminat ion of the galaxy clustering 
bias is at 10% level for low redshift galaxies [90L llOTt ll0S| . This may or may not be 
achievable at higher redshift, but it is not critical as BAO and WL, along with CMB, 
can determine the galaxy clustering bias quite well. 

We infer from the systematic errors (including extinction, photometry calibration, 
and seeing) of the SDSS angular galaxy power spectrum [§2 that 10~ 8 < N& ys < a few 
10~ 7 . For WL, the specified systematic shear noise is somewhat arbitrary and perhaps 
conservative, because its amplitude is relatively high compared to that of the shear 
power spectra. This partially compensates the fact that we restrict N^ ys to a very 
simple additive form. We also note that the effect of the simple noise N^ ys saturates 
quickly at N^ yg > a few 10~ 8 j^. In what follows, we assume the optimistic case by 
default, and the pessimistic results are presented in table 0] 

4. Forecasts 

The puzzle of dark energy has inspired a number of ambitious photometric and 
spectroscopic surveys that are proposed to unravel the mystery via multiple techniques. 
Among them, we choose to forecast cosmological constraints of BAO and WL for four 
multiband imaging surveys: DES, LSST, PS4, and SNAP in conjunction with CMB 
measurements by Planck, which is always included. It should be noted that all the 
four surveys propose to utilise one or two other techniques, SN and cluster counting, 
to constrain cosmology, and that BAO and/or WL may not even be their primary 
dark energy probe. 

We use the error product (EP) 

a{w p ) x a(w a ) = [a 2 (w )a 2 (w a ) ~ Cov 2 (w ,w a )} 1/2 , (15) 

where Cov(u>o,w a ) is the covariance between wq and w a , to assess dark energy 
constraints throughout this section. The minimum error of the dark energy EOS 
cr(wp) is achieved at the pivot expansion factor 

a p = 1 + Cov(w , w a )/o- 2 (w a ) (16) 

and is equal to the error on wq with w a held fixed [13, Il09) . This product 
is proportional to the area of the error ellipse in the WQ-w a plane and inversely 
proportional to the figure of merit used by the Dark Energy Task Force. 
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4.1. Surveys 

We model the underlying galaxy redshift distribution with 

ncxzV 2/z \ (17) 

which is in good agreemen t wit h photometry simulations based on degraded Hubble 
Deep Field North images 110]. The normalisation is adjusted to match the total 
surface density of galaxies n, and the galaxy distribution peaks at 2z*. Specifications 
of the four surveys are listed in table HJ 

Following [19j, we adopt a redshift-dependent rms shear 

7 rms = 0.18 + 0.042 z (18) 

for all the surveys. Smaller rms shear may be achieved with a subsample of well 
measured galaxies, and one can optimise the trade between 7 rms and rii to reduce the 
shape noise ^^n^ 1 f° r WL. It can also be beneficial to weight galaxies differently 
based on their measurement quality such as that of the shape noise and photo-z to 
optimise cosmological constraints. 



4-2. Dark energy constraints from BAO and WL separately 

The binning of the data affects the amount of information one can extract from the 
survey. Due to its broad kernel, WL only requires a small number of tomographic 
bins to fully reach its potential. For DES, five bins are sufficient jl^. We use 10 
equal-width WL bins from z p = to 3.5 for all the surveys. BAO, on the other hand, 
has narrow kernels that are defined by the galaxy distribution m, so that one can 
continue to extract information with finer bins until shot noise overwhelms the signal 
or the bin widths are less than the rms photo-z error. 

FigureElpresents the EP as a function of the number of BAO bins for each survey 
without the systematic noise, which can only reduce the number of bins required. We 
set the width of the bins to be proportional to (1 + z p ) to match the rms photo- 2: error. 
As expected, one can improve the dark energy constraint by using more tomographic 
bins until it saturates. The point of saturation depends on the depth of the survey 
and a z . For LSST (solid line) and SNAP (dotted line), the EP of 40 bins is only 7% 
worse than that of 90 bins, so, hereafter, we use 40 BAO bins for all the surveys. With 
uniform BAO bins (open circles) , the EP is more sensitive to the number of bins used. 
This suggests that it may be useful to optimise the tomographic binning for BAO. 
One possibility is to evenly distribute the (auto power spectrum) signal to noise ratio 
for each bin. 

BAO and WL, individually, are subject to different systematics and parameter 
uncertainties, but their combination can be substantially more robust. We first 
quantify in figure |3] the effect of the photo-z prior and fiducial value of the additive 
systematic noise on LSST BAO (left panel) and WL (right panel) using the EP. The 
BAO EP varies slowly with ap(5z) and Nj= ys , because LSST BAO is capable of self- 
calibrating the photo-z parameters to 10 _3 (1 + z) level (see figure and the BAO 
systematic noise is fairly low compared to the galaxy power spectra. WL is unable 
to constrain the photo-z parameters, so that it is sensitive to the photo-z prior on 
ap(Sz). The degradation to WL EP caused by the shear systematic noise saturates 
when N^ ys > a few 10 -8 , which is consistent with j^. Figure 03 demonstrates that 
to take the advantage of the WL technique one needs to know the photo-z error 
distribution accurately and have adequate control over the shear systematics. 
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Figure 2. The EP, u(w p ) X a(w a ), of BAO as a function of the number of 
tomographic bins for the four proposed surveys: DES (dash-dotted line), LSST 
(solid line), PS4 (dashed line), and SNAP (dotted line). The results are calculated 
without varying the number of parameters for the galaxy bias, photo-z bias, and 
photo-z rms, because they are defined in true-redshift space regardless the binning 
of galaxies in photo-z space. The bin widths are proportional to (1 + z p ) except 
that open circles are obtained with equal-size bins for LSST. The constraints are 
evaluated for the optimistic case in table but without the systematic error, 
which slightly reduces the number of bins needed to fully extract cosmological 
information from galaxy power spectra. 
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Figure 3. Left panel: Contours of the EP, ct(ui p ) X cr(w a ), of LSST BAO as a 
function of the prior of the photo-z parameters and fiducial value of the simple 
additive systematic noise. To reduce the dimension of the parameter space, we set 
the prior on the rms photo-z errors to 0"p(<r z ) = \/2o"p(<5z). The solid dot locates 
at the coordinates specified by the optimistic values of ap(8z) and Ni ys . Right 
panel: Same as the left panel, but for LSST WL. Notice that the EP of LSST 
WL varies by a factor of 5.9 across the ranges of the photo-z prior and systematic 
noise shown (from the lower left corner to the upper right corner), whereas the 
EP of LSST BAO changes by only 18%. If 7Vf ys = Af s 7 ys = and ct p (<5z) remains 
the same as that in the optimistic case, one could achieve EP = 0.0231 and 0.0068 
with BAO and WL, respectively. 



Cosmic tomographies 





1 


1 1 1 1 1 1 1 I 1 1 1 , 1 ] , 1 1 


0.1 




































b 




" .'/ • BA0 






\ '•.....•/ ° BAO Auto 


0.01 




\ / BAO + WL - 




"1 


\ / — - prior ; 
> ■ ■ 1 .... 1 .... 1 ... ~ 





1 2 




z 




z 



Figure 4. Left panel: Constraints on the galaxy clustering bias parameters from 
LSST BAO (with both the auto and cross power spectra, solid circles), BAO with 
the auto power spectra only (open circles), and joint BAO and WL (solid line). 
The 15% prior on each parameter is drawn in a thin dashed line. The constraints 
are tightest near z = 1, roughly corresponding to the peak of the overall galaxy 
distribution. Middle panel: Constraints on the photo-z bias parameters. The 
legends are the same as in the left panel with the addition of LSST WL in solid 
squares. The prior is taken to be ap(8z) = 0.05<r z = 0.0025(1 + z). Right panel: 
Same as the middle panel, but for the photo-z rms parameters. The prior is a 
factor of \/2 weaker than that on 8z. The galaxy cross power spectra are indeed 
crucial for self-calibrating the photo-z parameters and constraining the galaxy 
clustering bias. The oscillatory behaviour of WL constraints is caused by the 
offsets between the photo-z parameters and each tomographic bin. 

For BAO, the linear galaxy clustering bias b is degenerate with the linear growth 
function 67. Therefore, one cannot extract much useful information from the growth of 
the large-scale structure with BAO. It may seem contradictory from the left panel of 
figure^lthat the galaxy bias b (and, hence, the growth 67) can actually be determined 
to several percent with BAO (solid circles). This is made possible by normalising 
the matter power spectrum to the precisely measured CMB fluctuations so that the 
difference between the galaxy and extrapolated matter power spectrum at low redshift 
gives the galaxy bias directly. Even though the error of the galaxy bias can be as small 
as a few percent, it is still not a match to that of the distance from BAO [sil Eij . 
Without the galaxy cross power spectra (open circles) , the BAO results become several 
times worse, showing that the cross power spectra are useful for constraining the galaxy 
bias. Since WL is not affected by the galaxy bias, its combination with BAO (solid 
line) improves a(b) significantly. The tightest constraint on b occurs at a redshift that 
is lower than the peak of the overall galaxy distribution, because WL kernel peaks 
below the redshift of source galaxies. 

BAO and WL constraints on the photo-z parameters Sz and a z are shown, 
respectively, in the middle and right panels of figure 01 One sees that WL (solid 
squares) does not have much control on the photo-z parameters because of its 
broad kernel, while BAO, with both the auto and cross power spectra, is capable 
of calibrating the parameters to 10 _3 (1 + z) level. The galaxy cross power spectra 
play an important role in the BAO self-calibration of the photo-z error distribution; 
without them, the BAO constraints on the photo-z parameters become much weaker. 
Figure 21 thus explains the markedly different response to the photo-z bias prior 
between the BAO and WL EPs in figure It should be emphasised that the BAO 
self-calibration of the photo-z parameters cannot replace spectroscopic calibrations, 
because, without knowing how to faithfully parametrize the photo-z error distribution, 
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Figure 5. Complementarity between BAO and WL. The la wo-w a contours 
are shown for LSST BAO (dotted line), LSST WL (solid line), joint BAO and 
WL (shaded area), and the case of BAO plus growth information (dashed line), 
for which we fix the galaxy clustering bias parameters to their fiducial values. 
The last case is unrealistic, but it demonstrates the importance of the growth 
information. 

the self-calibration will be less useful. 

The photo- z constraints are much tighter than those of the galaxy clustering 
bias. The reasons are that the galaxy power spectra are more sensitive to the 
photo- z parameters than to the galaxy bias and that the dependence of the galaxy 
power spectra on the photo- z parameters is very unique. For instance, the derivative 
dlnnifijPf^ / db is roughly unity at z ~ 1 for all i and j, whereas AhinifijPf^ /dSz 
varies from —4 to 4 and behaves distinctively different from one pair of i and j to 
another (figure 

4-3. Joint analysis of BAO and WL 

Based on figure |3] and figure ^ one expects a great improvement of errors on 
cosmological parameters when BAO and WL are analysed jointly, which is indeed 
demonstrated in figure [3] The combination (shaded area) reduces the EP by more 
than a factor of 6 compared to that of LSST WL alone (solid line). This is more 
than the factor of ~ s/2 reported in |33j |. which does not include systematic errors 
or uncertainties in the photo- z error distribution. LSST BAO (dotted line) does not 
seem competitive to WL at all, if assessed by a(wo) and a(w a ) only. Since the errors 
in a(wo) and a(w a ) are highly correlated, the EP of BAO is merely a factor of 1.47 
larger than that of WL (see figure El and table 0J. The important issue here is not 
that which technique is better, but rather that one must analyse BAO and WL data 
(and possibly others) jointly to overcome the weakness of each. 

To understand the role of the growth function, we artificially fix the galaxy 
clustering bias parameters for BAO, which enables BAO to utilise the growth 
information. Since the growth function and galaxy bias are degenerate, the error 
of the growth function, if similarly parametrized, is the same as that of the galaxy 
bias. However, when the galaxy bias is fixed, the error of the growth function will 
be much smaller than that of the galaxy bias in figure 0] The resulting constraints 
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Figure 6. Left panel: The EP as a function of the rms photo-z error cr z 
for LSST BAO (dotted line), LSST WL (dashed line), and the two combined 
(solid line). We have applied the same tomographic binning regardless the value 
of cr z . With finer bins, BAO results can be further improved at small <r z , 
although WL and the combined results will remain less sensitive to a z . Right 
panel: The EP as a function of the normalisation of the primordial curvature 
fluctuation H 2 R at k = 0.05 Mpc" 1 . As A R increases, the galaxy number density 
fluctuation and weak lensing shear signal increase. If the signal is dominant 
over the noise (in a single mode), then the constraints on wo and w a will not 
change with the normalisation, whereas in the opposite case cr(u>o) and o(w a ) 
will be inversely proportional to the normalisation. To isolate the effect of the 
normalisation, we use the same maximum multipole of BAO, which is determined 
for A|j = 2.6 X 10 — 9 , for all normalisation cases. This means that the constraints 
with A 2 R < 2.6 x 10~ 9 are calculated from fewer multipoles than what would 
otherwise be included with the criteria A^(fc; z) < 0.4 and 40 < I < 3000. 



(dashed line in figure |5J) on wq and w a are, not surprisingly, as tight as those of BAO 
and WL combined, because the narrow kernel of BAO makes it possible to sample the 
cosmic density field at very fine intervals in redshift [see also|2l|. We note, however, 
that the determination of the (product of the galaxy clustering bias and) growth and 
distance are correlated. If the distance is poorly measured, the growth is not likely 
to be known very well either. Conversely, an accurate determination of the growth 
can improve the distance measurements, because the amplitude of the BAO and WL 
power spectra depend on the distance (not just distance ratios). 

Error forecasts inevitably depend on the assumed fiducial model. For example, 
because of its narrow kernel, BAO is sensitive to the rms photo-z error a Zl which 
can vary greatly with the galaxy population, filter design, survey depth, and so on. 
Another parameter of interest is the normalisation of the matter power spectrum, 
which is proportional to the normalisation A|j of the primordial curvature power 
spectrum at k = 0.05 Mpc -1 . Accompanying a lower estimate of the electron optical 
depth to the last scattering surface r = 0.09 and a slight tilt n s = 0.951, the recent 
analysis of the three-year WMAP date lowers the value of A|j by 25%. This reduces 
the amplitude of BAO and WL power spectra and can result in a weakening of the 
parameter constraints by up to 25% if the shot (shape) noise and systematic noise are 
dominant over the (auto power spectra) signal in a single mode. 

Figure quantifies the dependence of the EP on a z (left panel) and A^, (right 
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Figure 7. Left panel: Joint BAO and WL constraints on the dark energy EOS 
parameters wq and w a for DES (dash-dotted line), LSST (solid line), PS4 (dashed 
line), and SNAP (dotted line). Right Panel: Same as the left panel, but for the 
matter density LUm and curvature term Ck- Note that Planck alone constrains 
ujyn to ~ 1% (see table [3) and that SN and cluster counting can further improve 
the constraints. 



panel). The sensitivity of BAO (dotted line) and, to a lesser degree, WL (dashed 
line) EPs to a z is clearly seen, and the combination (solid line) is less sensitive to 
<j z . From the right panel, one sees that the BAO EP does not change much with the 
normalisation A|j. This means that it is not limited by noise, but rather by the photo- 
z error: with 40 bins whose widths are proportional to (1 + z p ), the bins are already 
narrower than the photo- z rms. The WL EP improves by 28%, i.e., roughly 14% each 
on w p and w a , as the normalisation increases from A 2 R = 2 x 1CP 9 to 2.6 x 10 -9 , 
indicating that the noise is more important for WL but not yet completely dominant. 



4-4- Forecasts 

We present wo~w a and w m -J7K If error contours of joint BAO and WL analysis for the 
four surveys in figure[7]and list the 1 a error bounds on all the cosmological parameters 
in table (optimistic case) and table ^ (pessimistic case) . 

Figure demonstrates that a deep and wide survey like LSST can tighten the 
constraints on dark energy and other cosmological parameters substantially. In its 
wide and shallow survey mode, PS4 can perform equally well as SNAP. Of course, the 
SNAP SN program will provide additional constraints that are fairly complementary 
to BAO and WL [e.g.,|2J. Likewise, the three ground-based surveys all plan to carry 
out SN and/or cluster counting to some degree. 

With w = —1, the combination of WMAP and the Supernova Legacy Survey 
data determines fix at a few percent level . The right panel of figure [7| shows 
that LSST BAO+WL and Planck can further reduce the uncertainty to 10~ 3 with 
our parametrization of the dark energy EOS, which may lead to interesting tests of 
inflation 0, 111 l| . It is also possible to measure the mean curvature from BAO and 
WL without assuming the dynamics and content of the universe [39|. Though less 
constraining, such a model independent method can provide a very useful check. 

Our forecasts for Planck are given in table [3] for reference. CMB is affected by 
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Table 3. Error forecasts with optimistic BAO and WL survey parameters 
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a The error of the dark energy EOS at the pivot epoch [equation 1151 1. a(w p ), is uncorrelated 
with that of w a . 



b When combining CMB with BAO and WL, we do not include the C MB c onstraints on wo 
and w a , which are sensitive to the late-time ISW effect on large scales |l!2| . 
c If the polarisation data at 2 < £ < 10 are included, Planck alone can achieve o-(r) = 0.0049 
and a(lnA 2 R ) = 0.010. 

d The running a = d In n B (k)/d In k at kf — 0.05 Mpc" 1 , and the primordial curvature power 
spectrum A^(fc) = A|(fcf)(fc/fc f ) na(fef) - 1+0 ' 5aln( ' s/fcf) is accurate to second order. 
c Without the self-calibration of photo- 2 parameters from the galaxy cross power spectra, the 
dark energy constraints would degrade to cr(w p ) = 0.021, a(wo) = 0.094, and a(w a ) = 0.22. 
One could achieve a(w p ) = 0.012, a(w ) = 0.048, and a(w a ) = 0.13, if Ni ys = NJ yB = 0. 



Table 4. Error forecasts with pessimistic BAO and WL survey parameters 
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dark energy mostly through the late-time ISW effect on large angular scales |ll2j . If 
we limit the multipoles to I > 40 for Planck, then its EP increases to 11. Since the 
properties of the cosmic density field on very large scales have not been measured 
very well yet, we neglect Planck constraints on the dark energy EOS parameters when 
combining it with BAO and WL. This treatment does not affect the combined results, 
especially those of all the three techniques, because wg and w a are better measured by 
BAO and WL. For non-dark-energy parameters (excluding Qk and a), the constraints 
are not so easily improved by BAO and WL except for LSST, which shows again that 
the deep photometry and wide area are not merely a technical challenge but rather 
a requirement for BAO and WL t o be competitive with future CMB measurements 
beyond the dark energy EOS |l!0| . 

As a result of its fine sampling of the density field in redshift space, BAO docs 
relatively well in measuring parameters (w m and a>b) that affect the features in the 
power spectrum. Since noises are more dominant for WL (see figure and figure EJ, 
BAO can perform better than WL for shallow surveys, especially in the pessimistic 
case (table BJ. BAO has another slight advantage that it probes the universe at a little 
higher redshift than WL does, because the lensing kernel peaks below the redshift of 
source galaxies. This advantage fades away when the survey is shallow and/or narrow, 
and it is a part of the reason that the relative strength of the BAO and WL constraints 
on Ok inverts from LSST to DES. For SNAP, its high galaxy number density (low 
shape noise) allows WL to remain significantly more competitive than BAO in the 
pessimistic case. 

Our errors of the running a for LSST WL and SNAP WL are considerably tighter 
than those in but the errors for DES BAO are larger than those forecasted for 
a /skv = 0.1 and z < 2.3 BAO survey in [40j. Since we do not include neutrino mass 
as does, and since assumes a survey with many more galaxies than DES, our 
results can still be consistent with theirs. Despite the differences, these forecasts all 
predict fairly tight constraints on a, which is expected to be of the same order as 
(n s — l) 2 from simple inflation models. With n s — 1 = 0.05 and measurements of a 
at 10 -3 level from future BAO and WL surveys, we will have another interesting test 
for inflation [Tij. 

Comparing the optimistic forecasts with pessimistic ones, we find that the joint 
BAO and WL results are fairly robust to the priors on non-cosmological parameters 
and to simple additive systematics. This clearly demonstrates the necessity of utilising 
multiple techniques jointly to extract cosmology information. 

5. Discussion and Conclusions 

We have shown that BAO and WL are highly complementary to each other especially 
in the presence of the uncertainties of the photo- z error distribution and inaccurate 
knowledge of the galaxy clustering bias. For non-Gaussian photo-z errors, more 
parameters will be needed to characterise the error distribution, which inevitably leads 
to weaker constraints on these parameters |86|. The impact of non-Gaussian photo-z 
errors on the forecasts of cosmological parameters will be investigated separately. 

We emphasise that one cannot rely solely on the galaxy power spectra for photo-z 
information. One may notice in figure 01 that the galaxy clustering bias and photo-z 
parameters are not constrained by the galaxy power spectra at both the low-redshift 
end, where very few power spectrum modes can be used, and the high-redshift end, 
where the data are poor. With long-tail photo-z error models, the unconstrained 
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errors at the two redshift ends can easily propagate to the whole redshift range and 
diminish the BAO ability to self-calibrate the photo-z parameters. In order to achieve 
the BAO self-calibration of photo-z errors, one must know how to parametrize the 
probability distribution of photo-z errors faithfully. In other words, it is crucial to 
map the photo-z error distribution and reduce its uncertainties as much as possible 
with massive spectroscopic calibrations and state-of-the-art photo-z algorithms. 

Even without knowing the true photo-z error distribution accurately, one can 
still use the tomographic BAO and WL power spectra to check for deviations from 
the assumed photo-z error distribution and to detect certain systematic errors. For 
example, a detection of the cross power spectrum between a very low redshift BAO 
bin and a high redshift one above the expected level could be a result of catastrophic 
photo-z errors between the bins or residuals of extinction correction. The effect of the 
former would be confined between the bins around the two redshifts, while that of the 
latter would be likely to spread over all bins. In any case, the cross power spectra 
arising from the two causes would have distinct shapes. Similarly, an unexpected 
detection of the cross power spectrum between a low redshift WL bin and a higher 
redshift BAO bin could be due to some systematics. 

We have assumed that the galaxy bias is not correlated with photo-z errors. In 
reality, different types of galaxies have different spectral features that affect their 
photo-z estimation and, meanwhile, do not cluster with the same strength. Therefore, 
photo-z errors are inevitably coupled with the galaxy bias through the galaxy type. 
More realistic forecasts should include multiple types of galaxies and assign separate 
galaxy bias and photo-z parameters for each type. Since the shot noise is relatively 
low for BAO, it is practical to split the galaxies into a few groups for better analyses. 
One may actually improve the photo-z constraints with multiple galaxy types |86| . 

One may worry about the break-down of the Limber approximation on large 
angular scales. The break-down of an approximation is not a loss of information, 
but, rather, a demand for proper means to extract the information. For parameter 
estimation, full calculations should replace the Limber approximation on large angular 
scales, but for error estimation, we can make more conservative cuts in £ to see the 
effect of completely ignoring the large-scale information. If we set the minimum 
multipole to £ m i n — 100, which removes the information of the broadband shape 
of the matter power spectrum (see figure the constraints on wq and w a degrade 
by less than 20% for LSST BAO and less than 13% for the joint analysis. In this 
test, BAO bins below z p ~ 0.23 have no multipole usable. BAO errors on wo and w a 
will double if one also excludes the first few acoustic peaks by setting £ m i n = 300. In 
contrast, WL errors will only increase by 25% with € m ; n = 300. These tests confirm 
that the constraining power of the broadband shape of the matter power spectrum is 
sub-dominant to that of the acoustic oscillations for the BAO technique 32] . 

On the other end of the scales, one may worry about our somewhat aggressive 
maximum multipoles £ meiK = 3000 for BAO and £ max — 2000 for WL. We find that 
a cut of £ max = 2000 for LSST BAO increases the error on wo and w a by less 
than 2%. This is not surprising, because the angular galaxy power spectra of the 
15 (z p > 1.7) BAO bins that are affected by this cut are practically featureless 
at 2000 < £ < 3000, i.e., these multipoles do not contribute much to distance 
measurements. With £ max = 2000 for BAO and £ max = 1000 for WL, the LSST 
BAO and WL joint constraints on wq and w a weaken by only 8%. Hence, we conclude 
that, given the survey parameters and assumptions about systematics, our forecasts 
are fairly robust against the range of multipoles included in the analysis. 
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Figure 8. Derivatives of the comoving angular diameter distance Da and linear 
growth function G with respect to uio, u> a , and Hk. When taking the derivatives, 
we hold the angular size of the sound horizon at the last scattering surface 8 B 
fixed. At z nj 0, the distance derivatives are effectively derivatives of the implicit 
Hubble constant with respect to the parameters. The linear growth function is 
normalised such that G(a) = a in an Einstein-de Sitter universe [as opposed to 
G(l) = 1 regardless the cosmological model]. 



There have been interesting discussions about dark energy constraints from 
measurements of the distance and the growth of the large-scale structure. The results, 
depending on how the distance and growth are separated, range from that they are 
equally powerful |64i Ir3^ | to that the constraints from the growth function are much 
weaker than those from the distance |2fJ, |62j • The relative strength of the distance 
and growth function in constraining dark energy depends on how sensitive the distance 
and growth function are to the parameters and how well they can be reconstructed 
from WL data. 

We plot the distance and growth derivatives with respect to Wq, w a , and Ok in 
figure [8] Since the growth function (open symbols) is less sensitive to the parameters 
than the distance (lines) and it is not so well determined as the distance with WL |20|. 
dark energy constraints from the growth are sub-dominant to those from the distance. 
However, this should not undermine the advantage of being able to measure the growth 
function (or, more appropriately, the amplitude of the matter power spectrum). To 
measure the distance with WL, one needs to know how the power spectrum amplitude 
changes with redshift, because the shear power spectrum depends on both the lensing 
geometry and growth function. In fact, the WL EP will increase by more than 10 
times, if one gives up the ability to measure the growth function by applying artificial 
bias parameters [with the fiducial value 6=1 and crp(ln6) = 0.2] to the shear power 
spectra [2lJ . Conversely, it is not possible to measure the growth function accurately 
if the distance is not, because one needs to know the scale k = £/Da accurately to 
estimate the amplitude of the power spectrum. 

The galaxy power spectra, being able to use the BAO features in the matter power 
spectrum, can better isolate the distance from (the product of the galaxy bias and) 
growth |44j, which partially explains the < 10% difference in the BAO constraints 
between the optimistic case [table [3] with ap(lnb) = 0.15] and pessimistic case [table U| 
with crp(ln6) = 0.3]. On the other hand, if the galaxy bias was known perfectly, then 



Cosmic tomographies 



22 




-1.2 -1 -0.8 

w o 

Figure 9. Error contours in the wo-w a plane with joint LSST BAO and WL 
analysis. Solid dots mark the fiducial values of wo and w a . All the four models 
have the same fiducial values of the other 9 cosmological parameters . The effects 
of wo and w a on the angular size of the sound horizon at the last scattering 
surface 9 a cancel each other along the thin dashed line. In other words, we require 
that the two implicit parameters f2x and h be the same for models on this line. 
The supcrgravity-inspired model at (wo,w a ) = (—0.8,0.6) has Qx = 0.61 and 
h = 0.57. 



one could determine the distance with both the amplitude and BAOs of the galaxy 
power spectra. This reduces errors in both the distance and growth function and leads 
to the tight dark energy constraints of the BAO+G case in figure [5J 

Figure |H1 also illustr ates why SN constraints on dark energy are sensitive to the 
flatness assumption The reason is that distance is much more sensitive to 

curvature then to the dark energy EOS. In addition, the (negative) distance derivative 
with respective to fix (dotted line) has a redshift dependence that somewhat resembles 
dlnD\/dw a (dashed line), so that a small freedom in Qk at the percent level can 
greatly degrade the constraint on w a . Similarly, the relatively tighter constraint on 
w compared to that on w a with all techniques can be explained by the relative 
amplitude of the derivatives. Since the distance derivative dlnDj^/dflK decreases 
more slowly with increasing redshift than dhi D^/dwo (solid line) and dlnZ?A/dw a , 
the degeneracy between dark energy and curvature can be lifted by extending the 
measurements to higher redshift. This explains the factor of ~ 3 difference in BAO 
(t(Ok) between the LSST survey and the shallower PS4 survey. The magnitude of the 
growth derivative with respect to fix (open circles) is at least 20 times larger than 
those with respect to wq (open squares) and w a (open triangles); the relative sign 
between the growth and distance derivatives with respect to fix is opposite to those 
with respect to wq and w a . Consequently, WL, being able to measure the growth 
function, can determine the curvature well even without high-redshift data, so that 
the PS4 WL constraint on fl^ is only moderately worse than that of LSST. 

Finally, we present the joint constraints of LSST BAO and WL for some 
interesting dark energy models in figure |5J All the models assume the same fiducial 
values for the other 9 cosmological parameters, i.e., essentially, the same distance to 
the last scattering surface. This requires modifications to the two implicit parameters: 
Six = 0.61 and h — 0.57 for the supergravity-inspired model at (wo, w a ) = (—0.8, 0.6). 
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For the models along the dashed line, fl x = 0.76 and h = 0.73 are preserved as well. 
In other words, w and w a are degenerate along the dashed line, if the distance to the 
last scattering surface is the only constraint on dark energy. Figure demonstrates 
that a joint analysis of LSST BAO and WL can indeed distinguish between some 
dark energy models that are quite degenerate to CMB or WL alone. Therefore, it 
is promising that, with multiple techniques as well as various statistics of the same 
galaxy and shear data, we will gain valuable insight to the mystery of dark energy 
from ambitious multiband imaging surveys. 
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Appendix A. Linear growth function 

Assuming that dark energy is spatially homogeneous, we can solve the the linear 
growth function G(a) from the linear perturbation equation 

G"(a)+^- {4 - a m (a) - n x (a) [1 + 3w(a)}} = B2&G{a), (A.l) 

where fl x (a) is the ratio of dark energy density to the critical density, and 
ft m (a) = a- 3 E- 2 {a) Q m , 

O x (a) = E- 2 (a)fl x exp 1^3 dlna' [l+w(a')] 

E 2 (a) = H^ 2 H 2 (a) = a~ 3 O m + a~ 2 (l - fl m - fi x ) + 

n x exp ^3 J dlna' [1 + w(a')] j , 

withfim = fi m (l) and fix = fix(l)- For a flat universe, (|A.1J| reduces to the from 
in |l!4| . With a matter-dominated universe at a -C 1, we set the initial condition to 
G(ai) — a l: where the starting expansion factor satisfies |G"(cti)| <C 1. The radiation 
component is dropped in the initial condition to follow the convention G{a) — a in an 
Einstein-de Sitter universe. This is appropriate since we are only concerned with the 
relative growth in matter era. The linear growth function may be normalised so that 
G(l) = 1 regardless the cosmological model in some applications. 



Appendix B. Nonlinear matter power spectrum 

Fitting formulae that are calibrated with A^-body simulations are widely 

used to map the linear matter power spectrum to nonlinear matter power spectrum. 
We use the PD96 formula 96] for convenience. However, modifications are necessary 
for the following reasons: (1) PD96 formula expects the logarithmic slope of the 
linear power spectrum to be monotonic, which is not satisfied with BAOs, and (2) it 
maps the linear wavenumber to nonlinear wavenumber, which would shift the scales 
of BAOs dramatically. To overcome these problems, we use a smooth fitting formula 
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|l!6j of the linear power spectrum without BAOs, A^ B L (fc), that otherwise matches 
the power spectrum with BAOs, A BAO L (fc), to calculate the nonlinear matter power 
spectrum, A^ B (fc). The final nonlinear matter power spectrum with BAOs is then 
A BAO (fc) = A| AO L (fc)A^ B (fc)/A^ B L (fc). We also note that one may use third-order 
pertu rbati on theory to calculate the nonlinear matter power spectrum with BAOs at 
z > 1 |ll7j ] , but there is no cle ar de monstration that the large error of the perturbative 
calculations at lower redshift |ll8| has been rectified. 
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